Skew scattering due to intrinsic spin-orbit coupling 
in a two-dimensional electron gas 
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We present the generalization of the two-dimensional quantum scattering formalism to systems 
with Rashba spin-orbit coupling. Using symmetry considerations, we show that the differential 
scattering cross section depends on the spin state of the incident electron, and skew scattering 
may arise even for central spin-independent scattering potentials. The skew scattering effect is 
demonstrated by exact results of a simple hard wall impurity model. The magnitude of the effect 
for short-range impurities is estimated using the first Born approximation. The exact formalism we 
present can serve as a foundation for further theoretical investigations. 



PACS numbers: 03.65.Nk, 71.70.Ej, 72.25.Rb, 73.50.Bk 
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I. INTRODUCTION 

When an electron scatters on an impurity atom with 
strong spin-orbit interaction (SOI) in a solid, the fea- 
tures of the process may depend on the spin state of 
the electron. The presence of the SOI may result in 
the asymmetry of the differential scattering cross sec- 
tion (commonly referred to as the skew scattering), even 
though the overall spin-dependent potential of the atom 
is central^* 2 In ferromagnetic metals, this particular scat- 
tering mechanism and the finite equilibrium polarization 
of the conduction electrons can lead to a finite poten- 
tial drop transverse to an applied electric field, even in 
the absence of magnetic field (anomalous Hall effect) 
Similar conditions in non-magnetic materials can induce 
a spin accumulation at the edges of the sample (extrinsic 
spin Hall effect) 

It is well known that the SOI can also be impor- 
tant in clean semiconductor samples lacking inversion 
sy mmet ryj 1 1 1 1 ' 1 2 For example, a finite electric field per- 
pendicular to the plane of a two-dimensional electron gas 
in a quantum well causes the spin splitting of the con- 
duction subbandst^ A widely used model Hamiltonian 
to describe this intrinsic spin splitting was proposed by 
Rashbaji 2 - 
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valid in the one-band effective-mass approximation. Here 
a x and a y denote the Pauli matrices, and the parameter 
a describes the strength of the spin-orbit coupling. The 
second term in (TJj) is usually referred to as the Rashba 
SOI or the Rashba term. 

Several studies have been devoted to the physical con- 
sequences of the interplay of the Rashba SOI and impu- 
rity scattering i 9 ' 13 i 14 i 15 i 16 i 17 i 18 The formation of a spin- 
polarized electron cloud around an isolated impurity has 
been predicted, resulting from the concurrent presence of 
the Rashba SOI and a constant in-plane electric fields 
The effect of the Rashba term on transport (intrinsic spin 



Hall effect) and polarization phenomena in disordered 
systems has been investigated using semiclassica l 9 ' 14 ! 15 
and quantu m 16 ! 17 ' 18 methods. In these studies, the role 
of impurity scattering in the presence of Rashba SOI is 
treated using various approximate models. However, it is 
necessary to go beyond these models if we want to reveal 
the details of the electron scattering characteristic of the 
presence of the Rashba coupling. 

In this work, we provide the generalization of two- 
dimensional scattering theory for the Rashba Hamilto- 
nian Ho via the S matrix formalism. Using only sym- 
metry considerations, we show that the Rashba term 
can induce skew scattering even if the scattering po- 
tential is central and spin-independent (e.g. an impu- 
rity ion or atom with negligible SOI). We demonstrate 
the skew scattering effect on the exactly solvable hard 
wall impurity mode L 19 ' 20 We prove that the effect ap- 
pears in the first Born approximation, which is a ma- 
jor difference compared to the conventional skew scatter- 
ing mechanising^ Based on this result, we propose a 
modification of the frequently used isotropic or s-wave 21 
approximation, in order to take into account the effects 
caused by the intrinsic spin-orbit coupling. 



II. GENERAL FORMULATION 

Consider the scattering problem governed by the 
Hamiltonian H = Ho + V , where the only assumption 
for the scattering potential V is to be zero outside a cir- 
cle of radius R. The cylindrical wave eigenfunctions of 
the Hamiltonian H with energy E are^ 2 - 



TH {d \{k T r)e-^/ 2 
H ( : d f(k T r)e^/ 2 
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where k = y / 2m*E/h 2 + k 2 Q , k T = 
am*/h 2 , t G {±1} is the helicity quantum number, 19 
and j G JT is the total angular momentum quantum num- 
ber. Here JJ = {...,-§, -A, ±, §,... }. d G {1,2}, and 



fji 1 ' 2 ) refers to the Hankel function of the first and second 
kind, respectively^ (h^) is an outgoing (incoming) 
cylindrical wave^ The magnitude of the radial particle 
current carried by the state ^ is independent of j, d and 
t. The eigenfunctions of Hq are eigenf unctions of H in 
the region r > R. Therefore, the elastic scattering of a 
single incoming partial wave hj T on the potential V is 
described by the wave function 

fr' 

valid outside the circle of radius R. The sum is for j' E JJ 
and t' e {±1}- The S matrix depends on the actual form 
of the scattering potential. Since the partial waves carry 
the same amount of current, the S matrix is unitary. In 
this section, we derive the relation between the S matrix 
and the differential scattering cross section. 

The plane wave eigenfunctions of Hq with energy E 
and helicity t G {±}, propagating in a direction ifi are 

<l> T , Vi (r,v)=r )r (v i y Krc ° a ^-V i \ (4) 

where 

VriVi) = -j= ("^ ) • (5) 
The partial wave expansion of the plane wave ([4]) is^ 

*™ = \{f T Y. l3+l/2[h ^ + h ^ e ~ iM - (6) 

T 3 

Using this expansion, the principal asymptotic form of 
the Hankel functions, and equation ([3]), it can be shown 
that the total wave function describing the scattering of 
the plane wave <j) TtVi asymptotically far from the scatterer 
is 

=0T^ + &% (7) 

with 

#$(r,¥0 = ^=J2 eikT ' r ^'(<P)fr'A<P^i)- (8) 

t' 

Here 

fM<P,<Pi) = ^=^e''^f)^e^^f), (9) 

and FjiS* = sQf* — 6jtjS T r T . In the following we will 
refer to the 2x2 complex matrix / as the scattering 
amplitude matrix. 

A plane wave with a given energy and propagation di- 
rection is not necessarily a helicity eigenstate. Coherent 
superpositions of the two helicity-eigenstate plane waves 
are described by the wave function 

07,¥» = ^ ^r.gtTri (10) 
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FIG. 1: (color online) Demonstration of the spin dynamics of 
plane waves a) with and b) without definite helicity (Eq. HQj). 
a) 7 = (1,0), b) 7 = (l,i)/y/2. In both figures, the black 
arrow denotes the direction of propagation (tpi = 0), the red 
arrows denote the local polarization vector, and the circle 
denotes the origin, tto is the polarization vector in the origin. 

where 7 = (7+,7_) T <E S*i(C 2 ) (the unit circle in the 
usual C 2 Hilbert space). The spin dynamics of such a 
superposed plane wave is essentially the same as in a 
Datta-Das spin transistor^ We show the spin dynamics 
of two examples in Fig. [1] 

In order to consider the scattering of the plane waves 
without definite helicity in (|10p . we introduce the 2x2 
complex matrix 

D(r,<p) = (e lfc +^+M,e 4fc -VM) ■ (11) 

Note that for every r and tp, D is unitary. With this no- 
tation, the scattered part of the wave function describing 
the scattering process of the superposed plane wave (fT0|) 
is 

r * 

(12) 

The differential scattering cross section is 

a diff ^, l p i ; 1 )=r\^ i (r^)\ 2 . (13) 

Since D is unitary, the differential cross section does not 
depend on r, as it is expected: 

CT d iff(vWi;7) =7 t [/(^Vi)] t /(^Vi)7- ( 14 ) 

In order to get a more transparent formula, we take the 
expansion of the scattering amplitude matrix using the 
unit matrix ctq and the Pauli matrices <j\ , 02 , 03 : 

3 

f(<P> fi) = X! fi) a m- ( 15 ) 

m=0 

The differential scattering cross section in terms of these 
coefficients u is 

Odiff (<p, W, 7) = c (<^ V>i) + V (<A <Pi) ■ p (7), (16) 
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with 

c = £ M 2 , (17) 

m— 

v = 2 Re(u*u) -i(u x u*), (18) 
P(7) = 7V7. (19) 

Here a = (a- x ,ay,a z ), therefore P(7) is a three- 
dimensional unit vector for arbitrary 7 € Si(C 2 ). In 
equation (fTr?)) the differential scattering cross section cor- 
responding to the angle ip is written in terms of the S- 
matrix (hidden in c and v), the propagation direction 
of the incoming plane wave ifi and the complex vector 
7 characterizing the superposed, non-helicity-eigenstatc 
plane wave. The result (|16[) is completely general, valid 
for arbitrary finite-range scattering potentials. This for- 
mula will play a central role in the derivation of the main 
result of our paper. 

To substitute the rather formal quantity 7 with a real 
physical quantity in (| 16[) . we can use the polarization 
vector of the superposed plane wave in the origin, 

7r o (7,</>i) = P(07,*n( r = )) (20) 

instead. Note that the P function is defined in (fT9|) . It 
can be shown that the following relation holds between 
the formal quantity P(7) and the real physical quantity 
7ro(7, </?i) = 

(0 cos ifi sin ifi \ 
sin^i -cos^i ■ P( 7 ). (21) 
-10 J 

This relation can easily be checked for the examples 
shown in Fig. [TJ Equation (fT6|) together with (j2Tj) ex- 
presses the differential scattering cross section corre- 
sponding to the angle if as the function of the plane wave 
propagation direction </?j, the S'-matrix and the polariza- 
tion vector of the incoming plane wave in the origin ttq. 

III. SIMPLE SCATTERING POTENTIALS 

The derivation so far has been completely general. 
Now we restrict our analysis to special scattering po- 
tentials. We call the V scattering potential simple, if 
it preserves the three fundamental symmetries of the 
Rashba Hamiltonian Hq: time reversal (ia y C, where C 
is the complex conjugation), rotation around the z axis 
( J z — —ihd<p + ha z /2) and the combined symmetry of 
real space reflection and spin rotation (a y P x , where P x 
is the spatial reflection with respect to the x axis). It 
is clear that the spin-independent central potentials are 
simple. 

For simple scattering potentials, the three symmetry 
operations above are compatible with H, which will re- 
sult in a remarkable simplification of the S matrix and, 
therefore, of Cdiff • It can be shown that the consequence 



of the time reversal, rotational and combined symmetries, 
respectively: 



qU'j) 




(22) 
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_ T > re i(j-j')ng(jj',-j) _ 
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As a consequence of these relations, the scattering am- 
plitude matrix - and hence every quantity derived from 
that - depends only on the scattering angle 6 = ip — ifi. 
The explicit form of / for simple scattering potentials: 



/t,t(0) = 
/-r,r(0) = 




Note that the diagonal (off-diagonal) elements of the 
scattering amplitude matrix / are even (odd) functions 
of the scattering angle 9. 

By definition, skew scattering is absent in the process 
if the differential cross section Cdiff in (|16p is an even 
function of the scattering angle 6 for every 7 6 Si(C 2 ). 
Using the symmetry properties of the components of the 
scattering amplitude matrix /, one can show that c and 
i>3 are even functions of 0. On the other hand, v% and 
v 2 are odd. It means that the absence of skew scattering 
is not provided by the symmetry properties of the total 
Hamiltonian even if the scattering potential is simple. 
This is the main result of our paper. 

To be a bit more specific, we can say that if the in- 
coming plane wave has a definite helicity, i.e. 7 oc (1,0) 
or 7 a (0,1), then using (p]) we get P(7) |j (0,0,1), 
therefore equation (JTSJ) and the even character of c and 
V3 implies the absence of skew scattering. On the other 
hand, if the incoming plane wave is a finite superposition 
of the two helicity eigenstates (i.e. both components of 
7 are finite), then the skew scattering effect arises if v\ 
or vi are finite. We note that obviously our symmetry 
considerations are not capable to tell whether v\ or V2 is 
finite or not - having a specific V potential in hand, we 
have to solve the scattering problem and calculate the 
elements of the ^-matrix in order to learn the answer. 



IV. HARD WALL IMPURITY MODEL 

In order to demonstrate the predicted skew scattering 
effect, we present exact results for a hard wall potential: 

I 00 otherwise 

This potential V is simple. We refer to Rcfs.Qj^ andl20l 
for the derivation of the elements of the S matrix. 

We focus on the low-energy properties of the scatter- 
ing, i.e. kR <C 1, because this limit corresponds to the 
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most widely studied short-range impurity models i 16 ' 17 ! 18 
If we consider the scattering of electrons at the Fermi 
energy, then it is realistic to set the parameter k so /k be- 
tween zero and 0.1. 25 Exact results for such parameter 
values are shown in Fig. [2l where we have plotted the 
angular dependence of the key quantities c, V\, V2 and 
t>3 defined in eqs. (|17p and (|18p . Apparently, v\ and vi 
does not vanish for finite fc so , therefore skew scattering 
is present in the process indeed. Other features of the 
results are summarized as: (i) the quantity c is practi- 
cally independent of fc so ; (ii) for k so — we have v = 0, 
therefore the differential cross section (which is equal to 
c in this case) is spin-independent, symmetric in 9 and 
approximately a constant function of 9; (iii) for finite k so , 
the magnitude of v\ is much smaller than that of V2 and 
V3; the latter ones appear to have the same magnitude for 
a given fc so , and their magnitude seems to scale linearly 
with fc so . 

We present the exact differential cross sections calcu- 
lated using (TIrJ)) in Fig. [31 corresponding to the two ex- 
ample plane waves of Fig. [T] The a dig of the plane wave 
with definite helicity (a) is symmetric, but for the super- 
posed plane wave (b) the skew scattering effect is clearly 
visible even for the realistic value of k so /k = 0.1. 

In order to understand the features (i) - (iii) of the 
exact results for the key quantities c and v (shown in 
Fig- we calculate them in the first Born approxima- 
tion for the simplest potential modeling short-range im- 
purites: V(r) = Kd(r), where S is the Dirac-delta and k 
represents the strength of the potential. For the incom- 
ing plane wave in (I10p . the scattered wave within this 
approach isi 2 . 



(28) 



where is the retarded Green's function of the Rashba 
Hamiltonian Hq. The exact form of the Green's function 
in position representation is know m 20 ' 22 Exploiting the 
simple form of our potential, we find 



< c 2( r ) 



K G+(r,O)0 7 , w (O), 



(29) 



where G^(r, r') is the position matrix element of G^. 
The first step to derive the scattering amplitude matrix 
is taking the |r| — > oo limit of the actual form of G^(r, 0), 

~ (sc) 

and calculating V'-y.vi- After that, one can derive the 
components of / using (TT2")) . With some algebra one gets 
the following results with respect to c and v: 



c(9) 
MO) 
v 2 (0) 

v 3 (9) 



(m*K) 2 
2irh 4 k 
0, 



Co, 



sin(#)co 
cos(#)c - 



k S o 

X' 



(30) 
(31) 
(32) 

(33) 



The qualitative similarity between these results and the 
exact ones in Fig [5] is remarkable. Apparently, the re- 
sults of this simple short-range impurity model grasp 
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FIG. 2: (color online) The quantities c, v\, V2 and vs (in units 
of R) determining the differential scattering cross section (see 
eq. (|16[0 as functions of the scattering angle 9 = if — ipi for a 
central hard wall potential. kR = 0.04. 



all the features (i), (ii) and (iii) of the exact results 
for kR <C 1 listed before. For k so — we recover the 
isotropic, spin-independent differential scattering cross 
section (Jdiff(^) = Co, as it is expected. Including only 
one further parameter k so /k, our results (|30H33[) provide 
a generalization of the one-parameter isotropic model for 
scattering in the presence of Rashba SOI. 



V. SUMMARY 

The knowledge of the properties of electron scatter- 
ing can be the starting point for further theoretical pre- 
dictions of impurity-related solid state phenomena. The 
theory of effects related to single, isolated impurities, like 
the Friedel oscillation, Landauer's charge dipole 2 ^ and 
the spin cloud predicted by Mal'shukov and Chu^ can 
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FIG. 3: (color online) Differential scattering cross sections 
exactly calculated in the hard wall impurity model, corre- 
sponding to the incoming plane waves of Fig.[T] a) 7 = (1, 0), 
b) 7 = (l,i)/V2. kR = 0.04 and k so /k = 0.1. 

be based on the knowledge of the scattering amplitude 
matrix /. In Boltzmann transport theory, the transi- 
tion probabilities are needed as an input information in 
the evaluation of the collision integral^ Hence the ex- 
act formalism presented in this paper can serve as a firm 



foundation of further theoretical investigations. The two- 
parameter model for short-range impurities provides a 
convenient tool to replace the isotropic approximation in 
systems with significant Rashba SOI. 

In conclusion, we have generalized the formalism of 
two-dimensional elastic quantum scattering to systems 
with finite Rashba spin-orbit coupling. Based on sym- 
metry considerations, we have shown that the differen- 
tial scattering cross section becomes spin-dependent and 
can show the skew scattering effect even if the scatter- 
ing potential is central and spin- independent. We have 
demonstrated the skew scattering by exact results of the 
hard wall impurity model. We derived the differential 
cross section in the first Born approximation for a Dirac 
delta scattering potential, and found remarkable similar- 
ity between these approximation and the exact results 
for low scattering energies. Using the simple formulas 
gained from the Born approximation, we proposed a two 
parameter model to substitute the isotropic or s-wave 
model of short-range impurity scattering in the presence 
of Rashba coupling. 
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Commission Contract No. MRTN-CT-2003-504574. 
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